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1. Introduction, motivations and notations 

The use of the discrete Fourier transform (DFT) is quite spread in many fields of physical 
sciences and engineering as for instance in signal theory. This chapter deals with a quadratic 
extension of the DFT and its application to quantum information. 

From a very general point of view, the DFT can be defined as follows. Let us denote (xg, X\, 
. . ., x&-\) a collection of d complex numbers. The transformation 

x = (x Q ,x- l ,...,x i -i) h-> x = (%*!,... ,*d-l) (!) 

defined by 

1 d-1 

x a = -)= £ e i2nan/d x„, ct = 0,1,..., d-1 (2) 
will be referred to as the DFT of x. 

Equation ((2) can be transcribed in finite quantum mechanics. In that case, x is often replaced 
by an orthonormal basis { | n) : n = 0, 1, . . ., d — 1} of the Hilbert space C d (with an inner 
product noted ( | } in Dirac notations). The analog of l|2) reads 

-1 d-1 

|«> = ~r E e i27ml/d \n), oc = 0,1, . . .,d - 1 (3) 
yd n =o 

Equation @ makes it possible to pass from the orthonormal basis { \n) : n = 0, 1, . . . , d — 1} to 
another orthonormal basis { | a) : a. = 0, 1, . . . , d — 1} and vice versa since 

(n\n f ) = 5(n,n') <^> (&\a') = 5(ot,a!) (4) 

The transformation ((3) defines a quantum DFT. In the last twenty years, the notion of quantum 
DFT has received a co nsiderable atten tion in connection with finite quantum mechanics and 
quantum information \ Vourdasl |2004|) . 

As an interesting property, the two bases {\n) : n = 0, 1, . . . , d — 1} and { 1 5) : a = 0, 1, . . . ,d — 
1}, connected via a quantum DFT, constitute a couple of unbiased bases. Let us recall that two 
distinct orthonormal bases 



B a = {\aa) :«. = 0,1,..., d-1} 



(5) 



and 



B b = {|fe|8) :B = 0,1,... ,d-l} (6) 
of the space C d are said to be unbiased if and only if 

V« = 0,l,...,d-l, Vj8 = 0,l,...,d-1 : |(««|6j8)| = -= (7) 

V d 

The unbiasedness character of the bases {\n) : n = 0, 1, . . .,d — 1} and { |&) : a = 0, 1, . . . ,d — 
1} then follows from 

(n\a) = A= e i2ncin/d => |(n|&)| = -^= (8) 
which is evident from {3). 

The determination of sets of mutually unbiased bases (MUBs) in C d is of paramount impor- 
tance in the the ory of information and i n quantum mechanics. Such bases are u seful in classi- 
cal information iCalderbank et alll 19971) , quan tum information iCerf et al.l,l2~002l) as well as for 
the constructio n of discrete Winner f uncti ons dGibbons et alj,l2004l) , the solution of the mean 
King problem dEnelert & Aharonovl.l200lh a nd the understanding of the Feynman path inte- 
gral formalism iTolar & Chadzitaskosl, 20091). It is well- known that the number N MUB of MUBs 
in C d is such that 3 < N^,, B < d + 1 iDurt et al.l,l201dl) . Furthermore, the maximum number 
N M ub = d + 1 is reached wh en d is a pr ime number or a power of a prim e number (Calder- 
bank et al., Il997t llvanovid, Il981 ; I Wootters & Fields! 19891) . However, when d is not a prime 



number or more generally a power of a prime number, it is not known if the limiting value 
N MU b = d + 1 is attained. In this respect, in the case d = 6, in spite of an e normous number 
of works it was not poss i ble to find more th an three MUBs (see for example jBengtsson et all 
l2007l;lBrierley & Weigeril2009l;lGrass]ll2005l) ). 

The main aim of this chapter is to introduce and discuss a generalization of the DFTs defined 
by (|2) and in order to produce other couples of MUBs. The generalization will be achieved 
by introducing quadratic terms in the exponentials in and ((3) through the replacement 
of the linear term an by a quadratic term ^n 2 + nn + Z, with n and £ in R. The resulting 
generalized DFT will be referred to as a quadratic DFT 

The material presented in this chapter is organized in the following way. Section 2 is devoted 
to the study of those aspects of the representation theory of the group SU(2) in a nonstandard 
basis which are of relevance for the introduction of the quadratic DFT. The quadratic DFT 
is studied in section 3. Some applications of the quadratic DFT to quantum information are 
given in section 4. 

Most of the notations in this chapter are standard. Some specific notations shall be introduced 
when necessary. As usual, 8 a j, stands for the Kronecker delta symbol of a and b, i for the pure 
imaginary, z for the complex conjugate of the number z, A f for the adjoint of the operator 
A, and I for the identity operator. We use [A, B\q to denote the ^-commutator AB — qAB 
of the operators A and B; the commutator [A, B] +1 and anticommutator [A, B]_i are noted 
simply [A, B] and {A, B}, respectively, as is usual in quantum mechanics. Boldface letters 
are reserved for squared matrices (Ij is the rf-dimensional identity matrix). We employ a 
notation of type \ ip), or sometimes \ip), for a vector in an Hilbert space and we denote (<p |i/>) 
and <p } (t/i respectively the inner and outer products of the vectors | ip) and 1 0) . The symbols 
ffi and © refer to the addition and subtraction modulo d or 2j + 1 (with d = 2j + 1 = 2, 3, 4, . . . 



depending on the context) while the symbol ® serves to denote the tensor product of two 
vectors or of two spaces. Finally N, N* and Z are the sets of integers, strictly positive integers 
and relative integers; R and C the real and complex fields; and Z/dZ the ring of integers 
0, 1, . . . , d — 1 modulo d. 

2. A nonstandard approach to SLT (2) 
2.1 Quon algebra 

The idea of a quon takes its origin in the replacement of the commutation (sign — ) and anti- 
commutation (sign +) relations 

a_fl+±fl+fl_=l (9) 
of quantum mechanics by the relation 

— qa + a_ = f(N) (10) 

where q is a constant and f(N) an arbitrary fonction of a number operator N. The introduction 
of q and f(N) yields the possibility to replace the harmonic oscillator algebra by a deformed 
oscillator algebra. For f(N) = 1, the case q = —1 corresponds to fermion operators (de- 
scribing a fermionic oscillator) and the case q = +1 to boson operators (describing a bosonic 
oscillator). The other possibilities for q and f(N) = 1 correspond to quon operators. We shall 
be concerned here with a quon algebra or ^-deformed oscillator algebra for q a root of unity. 
Definition 1. The three linear operators a_, a + and N a such that 

[a-,a+] q = I, [N a ,a + ]=a+, [N a ,a-] = -a-, {a+f = («_)* = 0, (N„) + = N a (11) 
where 

q = exp(^y fceN\{0,l} (12) 

define a quon algebra or q-deformed oscillator algebra denoted Aq(a-,a + ,N a ) or simply Aq(a). The 
operators fl_ and a+ are referred to as quon operators. The operators a~, fl+ and N a are called annihi- 
lation, creation and number operators, respectively. 

Definition 1 differs from the one by Arik and Coon jArik & Coonl[l976h in the sense that we 
take q as a primitive fcth root of unity instead of < q < 1. In Eq. (12) . the value = is 
excluded since it would lead to a non-defined value of q. The case fc = 1 must be excluded 
too since it would yield trivial algebras with a_ = a + = 0. We observe that for k = 2 (i.e., 
for q = —1), the algebra A_i(fl) corresponds to the ordinary fermionic algebra and the quon 
operators coincide with the fermion operators. On the other hand, we note that in the limiting 
situation where k — > oo (i.e., for q = 1), the algebra Aj(a) is nothing but the ordinary bosonic 
algebra and the quon operators are boson operators. For k arbitrary, N„ is generally different 
from fl + a_; it is only for k = 2 and — > oo that N a = a+a.-. Note that the nilpotency relations 
(a + ) k = (fl_) = 0, with fc finite, are at the origin of fc-dimensional representations of Aq{a) 
(see section 2.2). 

For arbitrary fc, the quon operators a_ and a + are not connected via Hermitian conjugation. It 
is only for fc = 2 or fc — > oo that we may take a + = («_ ) + . In general (i.e., for fc ^ 2 or fc -fr oo), 
we have (fl±) + ^ «=p. Therefore, it is natural to consider the so-called fc-fermionic algebra Eq 
with the generators a_, a + , a^_ = (fl+) + , si = and N a jDaoud e t al., 1998). The defining 



relations for E 9 correspond to the ones of Aq(a-,a + ,N a ) and Aq(a\,a^,N a ) complemented 
by the relations 

a_a^ — q~^a^a- = 0, — tpn+a + = (13) 

Observe that for k = 2 or k — > oo, the latter relation corresponds to an identity. The operators 
a.-, a + , a\_ and a_ are called fc-fermion operators and we also use the terminology fc-fermions 
in analogy with fermions and bosons. They clearly interpolat e between fermions and bosons. 
In passing, let us mention that the fc-fermions introduced in {Daoud et al .. 1998) share some 
common properties with the parafermions of order k — 1 discussed in {Beckers & Debergh , 
19901 ; iDurancU. Il993t iKharel, Il993l; iKhshevich & Plyushchayl Il999t iRubakov & Spiridonov . 



1988). The fc-fermions can be used for constructing a fractional supersy mmetric algebra of 
order k (or parafermionic algebra of order k — 1). The reader may consult toaoud et allll998|) 
for a study of the fc-fermionic algebra L ? and its application to supersymmetry. 

2.2 Quon realization of su (2) 

Going back to quons, let us show how the Lie algebra sw(2) of the group SU(2) can be gener- 
ated from two quon algebras. We start with two commuting quon algebras Aq (a) with a = x, y 
corresponding to the same value of the deformation parameter q. Their generators satisfy 
Eqs. frO and fl2t with a = x,y and [X, Y] = for any X in Aq(x) and any Y in Aq{y). Then, 
let us look for Hilbertian representations of Aq(x) and Aq(y) on fc-dimensional Hilbert spaces 
T x and Ty spanned by the bases {|«i) : «i = 0,1, . . .,k — 1} and { |«2) : n 2 = 0, 1, . . . ,k — 1}, 
respectively. These two bases are supposed to be orthonormal, i.e., 

(n 1 |ni) = J(n 1 ,ni) / (n 2 \n' 2 ) = 5{n 2 ,n' 2 ) (14) 

We easily verify the following result. 
Proposition 1. The relations 



x+\n{) = K + l), x + \k-l)=0 
x-\n{) = [ni\ 



1), x_|0)=0 (15) 



and 



y+|n 2 ) = [n 2 + l],|n 2 + l), y+|*-l)=0 

y-\n 2 ) = \n 2 -l), y-|0)=0 (16) 

N y \n 2 ) = n 2 \n 2 ) 

define k-dimensional representations of Aq(x) and Aq(y), respectively. In ([25} and $16[ , we use the 
notation 

VneN* : [n]q = ^f-^- = l + q + ... + q"-\ [0}q = 1 (17) 
1 q 

which is familiar in q-deformations of algebraic structures. 

Definition 2. The cornerstone of the quonic approach to su(2) is to define the two linear operators 



h= JN X (Ny + 1), V ra = S X Sy 



(18) 



with 

s y = y-r fl(N - N " ,/2 + ^ /2 [^r(y + )' c - 1 (20) 

In {29} and we take 



a e Z/dZ, cp r = nik - \)r, r 6 R (21) 

and the q-deformed factorials are defined by 

VneN* : [n],!=[l],[2] ? ...[n] ?/ [0],! = 1 (22) 

Nofe ife?f the parameter a might be taken as real. We limit ourselves to a in Z/dZ in view of the 

applications to MUBs. 

The operators h and v ra act on the states 

\ni,n 2 ) = \n\) ® \n 2 ) (23) 

of the A: 2 -dimensional Fock space T x ® Ty. It is straightforward to verify that the action of v ra 
on T x ® Ty is governed by 

v ra \k-l,n 2 ) = e i f' /2 \0,n 2 -l), n 2 ^ 

Vra\n x ,n 2 ) = q"^\ ni + l,n 2 -l), n x ^k-\, n 2 £ (24) 
v ra \n lr Q) = e^ /2 \ni + l,k-l), n^k-1 

and 

v ra \k-l,0) =e"f'\0,k-l) (25) 
As a consequence, we can prove the identity 

(v„) k = frl (26) 
The action of h on .F* ® .Fy is much more simple. It is described by 



h\n 1 ,n 2 ) = \jn x {n 2 + \)\n x ,n 2 ) (27) 

which holds for n\, n 2 = 0, 1, . . . , k — 1. Finally the operator v ra is unitary and the operator h 
Hermitian on the space T x ® Ty. 

We are now in a position to introduce a realization of the generators of the non-deformed Lie 
algebra su(2) in terms of the operators v ra and h. As a preliminary step, let us adapt the trick 
used by Sc hwinger in his a pproach to angular momentum via a coupled pair of harmonic 
oscillators iSchwingen. 11965ft . This can be done by introducing two new quantum numbers / 
and M 

J = I ( Wl + n 2 ) , M = - (m - nz) (28) 



and the state vectors 

\},M) = \n lr n 2 ) = \J + M,J-M) => (J,M\j',M') = 6jj,S MM , (29) 

Note that 

j=l(k-l) (30) 

is an admissible value for /. We may thus have j = 4, 1, |, . . . (since = 2, 3, 4, . . .)• For 
the value /' of /, the quantum number M can take the values m = j r j — 1, . . . , —j. Then, let us 
consider the (2/ + 1) -dimensional subspace e(j) of the fc 2 -dimensional space T x ® T-y spanned 
by the basis 

%+i = { I// m) : w = - 1, . . . , -/'} (31) 
with the orthonormality property 

{j,m\j,m') = &„ hm , (32) 

We guess that e(j) is a space of constant angular momentum As a matter of fact, we can 
check that e(j) is stable under ?/ and i> rB . 

Proposition 2. 77/e act/on o/ the operators h and v ra on e(j) is given by 

h\j,m) = J(j + m)(j — m + 1) \j,m) (33) 

v ra \j,m) = S m / 2 ^\j, -j) + (1 - V/)<7°'~ m)a |7> + 1) (34) 

tyfeere q is given by |12} with k = 2j + 1, r 6 E and a 6 Z/(2j + 1)Z. 
It is sometimes useful to use the Dirac notation by writing 

h = £ ^(J + m)(j-m + m)(/, m| (35) 

m=— / 

7-1 

v ra = e l2n ' r \j,-j){j,j\+ £ ^- m ) fl |;,m + l)(/,m| (36) 

m=-j 

(v m y = e - i2 ^\j,j)(j,-j\+ £ q -(i- m +V*\j,m-l)(j,m\ (37) 

m=-j+\ 

It is understood that the three preceding relations are valid as far as the operators h, v ra and 
(v ra ) f act on the space e(j). It is evident that h is an Hermitian operator and v ra a unitary 
operator on e{j). 

Definition 3. The link with sw(2) can be established by introducing the three linear operators j + , j- 
and j z through 

j+ = hv ra , ]- = (v ra f h, j z = j \ll 2 - (Vraf h 2 v rc ^ (38) 

For each couple (r, a) we have a triplet (]+,]-, j z )- It is clear that j + and /_ are connected via Hermitian 
conjugation and j z is Hermitian. 



Proposition 3. The action of j+, j- and j z on e(j) is given by the eigenvalue equation 

iz\i' m ) — m \i' m ) (39) 

and the ladder equations 

j + \j,m) = ,(;-«+s-i/2)y (j _ m )(j + m + i) \j t m + i) (40) 
j-\j,m) = q -U-»'+s+W)«^Q + m )(y _ m + X )\],m - 1) (41) 

where s = 1/2. 

For a = 0, Eqs. 1(39), ( f40t and fiTt give relations that are well-known in angular momentum 
theory. Indeed, the case a = corresponds to the usual Condon and Shortley phase conven- 
tion used in atomic and nuclear spectroscopy. As a corollary of Proposition 3, we have the 
following result. 

Corollary 1. The operators j + , /_ and j z satisfy the commutation relations 

[/*/+]=/+, [/*/'-] = -/'-, ]]+,]-]= 2} z (42) 
and thus span the Lie algebra of SU(2). 

The latter result does not depend on the parameters r and a. The writing of the ladder opera- 
tors /+ and /_ in terms of h and v ra constitutes a two-parameter polar decomposition of the Lie 
algebra sw(2). Thus, from two (/-deformed oscillator algebras we obtained a polar decompo- 
sition of the non-deformed Lie algebra of SU (2). This decomposition is an alternative to the 
polar decomposition s obtained independently in dChaichian & Ellinasl, ll99Ct iLevy-Leblondl, 
ll973l:IVourdaIll990|) . 

2.3 The {j 2 ,v ra } scheme 

Each vector \j, m) is a common eigenvector of the two commuting operators j z and 

j 2 = 2 0'+;'- + ;'-;'+) + ;1 = /+/'- + 73(73 - 1) = 7-7+ + 73(73 + 1) (43) 

which is known as the Casimir operator of su(2) in group theory or as the square of a gener- 
alized angular momentum in angular momentum theory. More precisely, we have the eigen- 
value equations 

j 2 \j,m) = j(j + l)\j,m), j z \j,m) = m\j,m), m =],]— 1, —j (44) 

which show that j and m can be interpreted as angular momentum quantum numbers (in 
units such that the rationalized Planck constant h is equal to 1). Of course, the set {j 2 ,j z } is a 
complete set of commuting operators. It is clear that the two operators j 2 and v ra commute. 
As a matter of fact, the set {j 2 ,v ra } provides an alternative to the set {j 2 ,j z } as indicated by 
the next result. 

Theorem 1. For fixed j (with 2j 6 N*), r (with reR) and a (with a 6 Z/(2; + 1)Z), the 2] + 1 
common eigenvectors of the operators j 2 and v ra can be taken in the form 



(45) 



ex P ( ) 



where 

The corresponding eigenvalues are given by 

j 2 \ja;ra) = j(j + l)\ja;ra), v ra \ja;ra) = q^ r+a ^ K \ia;ra) , a. = 0,1,. . .,2/ (47) 

so that the spectrum of v ra is nondegenerate and {j 2 ,v ra } does form a complete set of commuting 
operators. The inner product 

{joc;ra\jP;ra) = 6 Kr p (48) 

shows that 

Bra = {\w,ra) : a = 0,1,..., 2;} (49) 

is a nonstandard orthonormal basis for the irreducible matrix representation ofSU(2) associated with 
j. For fixed j, there exists a priori a (2j + l)-multiple infinity of orthonormal bases B ra since r can have 
any real value and a, which belongs to the ringl^l (2j + 1)Z, can take 2j + 1 values (a = 0, 1, . . . ,2j). 
Equation (45) defines a unitary transformation that allows to pass from the standard orthonor- 
mal basis B 2 y + i, quite well-known in angular momentum theory and group theory, to the 
nonstandard orthonormal basis B ra . For fixed r and a, the inverse transformation of (45) is 

1 2 i 

\j, m ) = q -(i+'n)(j->n+l)a/2+jmr^_ £ I? -(/+'«)«|y a;rfl ) ; m = j,j-l,..., -j (50) 

which looks like an inverse DFT up to phase factors. For r = a = 0, Eqs. (45} and (50) lead to 
|/«;00) = -j2= £ ?G +B >«|/,m), « = 0,1,...,2; (51) 

V l l + 1 m=-j 

« |;,m)=-==£ 9 -° +m) "l/«;00), m=;,;-l (52) 

V 2/ + 1 a=0 

Equations (ST) and (52) correspond (up to phase factors) to the DFT of the basis B 2 y + i and its 
inverse DFT, respectively. 

Note that the calculation of (/«; ra \jB; sb) is much more involved for (r ^= s, a = b), (r = s, a ^ 
b) and (r ^= s, a ^ b) than the one of (ja;ra\jf};ra) (the value of which is given by (48)). For 
example, the overlap between the bases B ra and B sa , of relevance for the case (r ^= s, a = b), is 
given by 

(icL-raUB-sa) = — Sln[/(S ~ r) + - ~ ^ (53) 

The cases (r = s, a b) and (r ^ s, a b) need the use of Gauss sums as we shall see below. 
The representation theory and the Wigner-Racah algebra of the group SU(2) can be developed 
in the {j 2 , v ra } quantization scheme. This leads to Clebsch-Gordan coefficients and (3 — jot) r a 



symbols with properties very different from the ones of the usual SU(2) D (1(1) Clebsch- 
Gordan coefficients and 3 — jm symbols corresponding to the {j 2 ,jz} quantization scheme. 
For more details, see Appendix which deals with the case r = a = 0. 

The nonstandard approach to the Wigner-Racah algebra of SU(2) and angular momentum 
theory in the {j 2 , v ra } scheme is especially useful in quantum chemistry for problems involv- 
ing cyclic symmetry. This is the case for a ring-shape molecule with 2/ + 1 atoms at the ver- 
tices of a reg ular polygon with 2/ + 1 sides or for a one-dimensional chain of 2j + 1 spins of 
i -value each jAlbouy & Kiblej,l2007l) . In this connection, we observe that the vectors of type 
\ja; ra) are specific symmetry-adapted vectors. Symmetry-adapted vectors are widely used in 
quantum chemistry, molecular physic s and condensed matt er physics as for instanc e in ro- 
vibra tional spectroscopy of molecules dChampion et aljjl977l) and ligand-field theory jKiblerl 
Il968h. However, the vectors \ia;ra) differ from the symmetry- adapted vectors considered in 
{Champion et alJ,ll977tlKiblerl,[l96a|Patera & Winternitzlll976l) in the sense that v ra is not an 
invariant under some finite subgroup (of crystallographic interest) of the orthogonal group 
0(3). This can be clarified as follows. 

Proposition 4. From ( f36t , it follows that the operator v ra is a pseudo-invariant under the cyclic 
group C2j+\, a subgroup of SO(3), whose elements are the Wigner operators Pr(^) associated with the 
rotations R{f), around the quantization axis Oz, with the angles 

p = Q,\,...,2j (54) 



r 2/ + l 

More precisely, v ra transforms as 



P R ( f) Vra [P R(f) ) =e-" p v ra (55) 
Thus, v ra belongs to the irreducible representation class ofC2j+\ of character vector 

X {2 ' ] = (hq- 1 ,...,q- 2 ') (56) 

In terms of vectors ofe(j), we have 

PR(<p)\wra)=qiP\jFj;ra), p = aQp (57) 

so that the set {\joc;ra) : a = 0, 1, . . . ,2j} is stable under Pr^- The latter set spans the regular 
representation ofC2j+\- 

2.4 Examples 

Example 1: The ] = \ case. The eigenvectors of v ra are 

\\a-,ra) = ^e in ^ /2 - r/i+K \l,l) + ^e inr,i \\,-l), a = 0,l (58) 
2 1 y^ '2 2' V2 2 2 7 x ' 

where r6E and a can take the values a = 0, 1. In the case r = 0, Eq. {58} gives the two bases 

1„ ™x 1 (A lv ,1 M ,1, ™x 1 /,1 1, ,1 1. 

. |-0; 00) = -j= I 

and 



3oo : |_0;00> = [\-,-) + \- r - 2 )), \~ 2 im = --^ \\-,-) - |-,- ? ) ) (59) 



Bor : l^;0i)^(|i,i)-4,-i)), \\^ = -^ 2 (\\,\)+i\\,-\)) 



The bases d59t and l l60t are, up to phase factors, familiar bases in quantum mechanics for 
|-spin systems. 

Example 2: The / = 1 case. The eigenvectors of v ra are 

\lu;ra) = -±=q r (q a+2a - 2r \l,l) + q e+el ~ r \l l Q) + a = 0,1,2 (61) 

where r 6 R and « can take the values a = 0, 1, 2. In the case r = 0, Eq. (|6T) gives the three 
bases 

B 00 : 1 10; 00) = (|1,-1) + |1,0> + |1,1» 

|11;00) = -^(|1,-1) + < 7 |1,0)+< 7 2 |1,1)) (62) 
1 12; 00) = ±=(\\,-\)+q 2 \l,0)+q\l,l)) 

% : |10;01) = ^=(|1,-1) + 9 |1,0)+*|U» 

1 11; 01) = -^(|1,-1)+^ 2 |1,0) + |1,1)) (63) 
1 12; 01) = (|1,-1) + |1,0)+^ 2 |1,1)) 

10; 02) = (|1/-1) +9 2 |1.0) +<7 2 |1,1)) 

1 11; 02) = -== (|1, — 1> + 1 1,0) +<j|l,l)) (64) 
v3 

1 12; 02) = -= (|1,-1) +q\l,0) + |1,1)) 
v3 

It is worth noting that the vectors of the basis Bqo exhibit all characters 

x («) = (l,f,q 2tt \ , a = 0,1,2 (65) 



: '02 



of the three vector representations of C3. On another hand, the bases Bqi and B02 are connected 
to projective representations of C3 because they are described by the pseudo-characters 

Xi ] =(Xq l+a ,q 1 -*), a = 0,1,2 (66) 

and 

xP=(W + W~"). « = 0,1,2 (67) 

respectively. 



3. Quadratic discrete Fourier transforms 

We discuss in this section two quadratic extensions of the DFT, namely, a quantum quadratic 
DFT that connects state vectors in a finite-dimensional Hilbert space, of relevance in quantum 
information, and a quadratic DFT that might be of interest in signal analysis. 

3.1 Quantum quadratic discrete Fourier transform 

Relations of section 2 concerning SU(2) can be transcribed in a form more adapted to the 
Fourier transformation formalism and to quantum information. In this respect, let us intro- 
duce the change of notations 

d = 2j + l, n = j + m, \n) = \j,-m) (68) 

and 

\aa;r) = \ja;ra) (69) 

so that (49) becomes 

B ra = {\aa;r) : a = 0,1,.. .,d- 1} (70) 

(Note that d coincides with the dimension k of the spaces T x and Ty of section 1.) Then from 
Eq. l |45t , we have 

\aot- r) = (? (rf-D 2 'V4 1 £ ^n(d-n)«/2+n[«-(d-l)r/2]| d _!_„), a = , 1, . . . , d - 1 (71) 
vd n=0 

or equivalently 

\a*;r) = ,(i-i)V4 1 £ 9 (d-l-»)(»+i)a/2+(i-l-n)[*-(d-i)r/2]| n)i a = 0/J d _ 1 (72) 

V d „=o 



where 

The inversion of (71) gives 



9 = exp ( ) 



|d-l-n) =(? -»(rf-")«/2-(rf-i) 2 >V4+«(rf-i) r /2 1 rf v- 1 i? -««| flfl ,. r ^ „ = 0,l,...,d-l (74) 

Vd x=0 

By introducing 

(F ra ) M = J_ q "(d-">n+(d-lfr/i+n[ a -(d-l)r/2] r fi,0i = 0,1, ... ,d — 1 (75) 

v d 

equations (71} and f74) can be rewritten as 

rf-i 

|«*;r) = £ (F ra ) na \d-l-n), oc = 0, 1, . . . ,d - 1 (76) 

11=0 



and 



d 1 

\d-l-n) = £ (F ra ) n J«a;r), n = 0, 1, . . . ,d - 1 (77) 

a=0 



respectively. For r = a = 0, Eqs. (76) and f77l yield 



-i d-l 

|0a;0) = -j= ^e i2nan/d \d-l-n), oc = 0,1, . . .,d - 1 



1 



rf-i 



^ |d-l-n) = ^ ^e- ,27I " a/rf |0a;0), » = 0, 1,. . .,d - 1 (78) 

yd a=0 

which corresponds (up to a change of notations) to the DFT described by @. For a ^ 0, 
Eq. (76) can be considered as a quadratic extension (quadratic in n) of the DFT of the basis 
{ |n) : n = 0, 1, . . . ,d — 1} and Eq. {77} thus appears as the corresponding inverse DFT. This 
can be summed up by the following definition. 
Definition 4. Let H ra be the d x d matrix defined by the matrix elements 

(H ra ), m = 1 1? (d-l-n)(n+l)«/2+(d-l)V/4+(d-l-»)[«-(d-l)r/2l ; „, a = ,1, . . . , i - 1 (79) 

v d 

where, for a fixed value of d (with d 6 N \ {0, I}), r and a may have values in ]R and "Z./d"Z., 
respectively. In compact form 

(H ra ), m = ^f mvld (80) 

with 

v = -i(d-l) 2 r+i(d-l)a + (d-l)a- ^[2a + 2a-da- {d-\)r]n- ^an 2 (81) 
The expansion 

d-l 

\aoc;r) = £ (H ra ) na \n), oc = 0,1,. . .,d - 1 (82) 

defines a quadratic quantum DFT of the orthonormal basis 

B d = {\n) :n = 0,\,...,d-l} (83) 

This transformation produces another orthonormal basis, namely, the basis B ra (see Eq. {70}). The 
inverse transformation 

d-l 



\n) = £ (H ra ) n Jfla;r), n = 0,1, . . .,d - 1 (84) 

ct=0 



gives back the basis B d . 



For fixed d, r and a, each of the d vectors \aoc; r), with a = 0, 1, . . .,d — 1, is a linear combination 
of the vectors 0), |1), . . . , \ d — 1). The vector \aa;r) is an eigenvector of the operator 



d-2 



v ra = e iK{ - d - x >\d - 1)(0| + £ q [d ~ X ~ n)a \d -2-n)(d-l-n\ 

ii=0 



or 



rf-1 

z, rfl = e M<i-l)r| d _ !) ( | + £ 9 ™| B - l)(n\ 

n=\ 

(cf. Eq. (36)). The operator v ra can be developed as 

Vra = e in ^ d - 1 >\d-\){0\ +<? fl |0)(l| +<? 2fl |l)(2| +... + q^- 1 >\d-2){d-\\ 
Then, the action of v ra on the state \n) is described by 

v ra \n) = 8 nM e in ^ r \d - 1) + (1 - S nt0 )c," a \n - 1) 
(cf. Eq. (34)). Its eigenvalues are given by 

v m \acc;r) = q^-^^'^a^r), a = 0,1,. . .,d - 1 

(cf. Eq. (47)). 



(85) 



(86) 



(87) 



(88) 



(89) 



3.1.1 Diagonalization of v, •,, 

Let V ra be the d x d unitary matrix that represents the linear operator D rfl (given by (87) ) on 
the basis B d . Explicitly, we have 



/ <? fl 

q 2a 



^ e i7t(d-l)r o o 



\ 


^(rf-i)a 
o J 



(90) 



where the lines and columns are arranged in the order 0, 1, . . . , d — 1. Note that the nonzero 
matrix elements of Vqu are given by the irreducible character vector 



x (") = (l, i? fl ,.--,<? (d ~ 1)fl ) 



(91) 



of the cyclic group Q. 

Proposition 5. The matrix H ra reduces the endomorphism associated with the operator v ra . In other 
words 



(H ra ) + V ra H ra = (? ( d - 1 )('-+«)/ 2 
in agreement with Eq. (47). 



fq° o 

q- 1 

Vo 



5 ] 

-(d-i)J 



(92) 



Concerning the matrices in {90) and {92}, it is important to note the following convention. 
According to the tradition in quantum mechanics and quantum information, all the matrices 
in this chapter are set up with their lines and columns ordered from left to right and from 
top to bottom in the range 0, 1, . . . ,d — 1. Different conventions were used in some previous 
works by the author. However, the results previously obtained are equivalent to those of this 
chapter. 

The eigenvectors of the matrix V ra are 



d-1 



(p(aa; r) = £ (H ra )„« <pn, * = 0, 1, . . . ,d - 1 

n=0 

where the <p„ with n = 0, 1, . . . , d — 1 are the column vectors 



(93) 



/ 2 \ 




l 



<pd-l 







w 



(94) 



representing the state vectors 0), 1 1), . . . , \ d — 1), respectively. These eigenvectors are the col- 
umn vectors of the matrix H ra . They satisfy the eigenvalue equation (cf.|89t 



V ra </>(aa;r) = 9 (d_1)(r+fl)/2 -*0(«a;r) 



(95) 



with ct = 0, 1, 



3.1.2 Examples 

Example 3: The d = 2 case. For d = 2, there are two families of bases B ra : the B,-q family and 
the B r i family (a can take the values a = and a = 1). In terms of matrices, we have 



H r 



(f 



«/2-r/4 



V2\ f 



r/A 



? fl/2-r/4 



q" 
q r 



The matrix V ra has the eigenvectors (corresponding to the basis £>, «) 

<p(aa;r) = -^(< 7 fl/2 - r/4 +> + q r/A <Pi), « = 0,1 
v2 



where 



For r = 0, we have 



(96) 



(97) 



(98) 



1 

1 



the eigenvectors of which are (cf. {97}) 
^(00;0) = -L(^+^) = ^(} 



^(01;0) = i=^-^) = -^(_ 1 1 



(99) 



(100) 



and 



*(iO;0) = -^(fr + ito) = -^ry, 0(ll;O) = ^(0i-z0 o ) = --^Q (101) 

which correspond to the bases Boo an d Bqi, respectively. Note that l llOOt and l llOlt are, up to 
unimportant multiplicative phase factors, qudits used in quantum information. 
Example 4: The d = 3 case. For d = 3, we have three families of bases, that is to say B y q, B r \ 
and B r 2, since a can be 0, 1 and 2. In this case 

(q a - r q a+2 - r q a+l - r \ / q a \ 

H ra = -L I q" q a+l q a+1 \, V ra = ( q 2a ) , ? = e C7t/3 (102) 







q" 




), v.= | 







'o) 












0(aa; r) = -L^ (q a+2K - 2r cp + q a+a - r (Pi + 2 ) , « = 0, 1, 2 (103) 



and V ra admits the eigenvectors (corresponding to the basis B r « 
1 

where 

/1\ /0\ /0\ 



00 = ^0 J , 0j = ^1 j , 2 = ^0j (104) 

In the case r = 0, we get 

/0 1 0\ /0 q 0\ /0 (? 2 0\ 

Vbo = 1 , % = ? 2 , V 02 = q) (105) 
\1 0/ \1 0/ \1 0/ 

The eigenvectors of Vbo/ ^bi an d Vfe follow from Eq. j!03t . This yields 

0(00; 0) = -j= {<p2 + 01 + <po) 



0(01; 0) = -j= ( 2 + #1 + <r>o ) (106) 



0(02; 0) = (02 + r0i+#o 



or 

„2 



0(OO;O) = -L f 1 J , 0(01; 0) = ±= U J , 0(02; 0) = -L ( (107) 



corresponding to Bqo, 



0(10; 0) = i(0 2 + #i + #o) 

0(11; 0) = (fcj + ,2^ + ^ (108) 
1 

7! 



0(12; 0) = -J= U 2 + cpi+q 2 4 



or 



*(io;0) = ^ 

corresponding to Bqi, and 



<p(12; 0) 



</>(20;0) = 



tf>(22; 0) 



1 

1 



</>2 + (? 2 </>i + <7 2 </>o 



0(21; o) = — (^2 + ^ + ^) 



v5 



(<fe + q<pi + 4>o) 



l 

75 



(109) 



(110) 



,(20;0)^ ft 



, ( 21;0)^ J 



4>(22; 0) 



1 

V3 



(111) 



corresponding to Bq2- Note that 4107b , 4109t and l llllt are, up to unimportant multiplicative 
phase factors, qutrits used in quantum information. 



3.1.3 Decomposition of V ra 

The matrix V ra can be decomposed as 



V ra = P r XZ" 



where 





( l 





. 
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. 







Pr = 
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e in(d- 





and 
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• °^ 




f 1 
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• \ 










1 . 


. 
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X = 










, z = 








<? 2 . 
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\l 
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• o) 









. 





The matrices X and Z can be derived from particular V ra matrices since 

X = V 00 , Z = (V 00 ) f V 01 
which emphasize the important role played by the matrix V ra . 



(112) 



(113) 



(114) 



(115) 



The matrices P r , X and Z (and thus V ra ) are unitary. They satisfy 

V ra Z = qZVrz (116) 
V ra X = ^- fl XV ra (117) 

Equation l !116t can be iterated to give the useful relation 

(V ra )"'Z" = 9 "»Z B (V r ,) m (118) 

where m,n E Z/dZ. Furthermore, we have the trivial relations 

e -m(d-l) r{ y to) d = Z d = I d (119) 

More generally, we can show that 

Vn 6 Z/dZ : (V ra )« = q- ni ~ n -^ a/2 (\ tQ ) n Z an (120) 

Consequently 

(V ra ) d = e in(d-l)(r+a) ld (m) 
in agreement with the obtained eigenvalues for V ra (see Eq. (95)). 
3.1.4 Weyl pairs 

The relations in sections 3.1.1 and 3.1.3 can be particularized in the case r = a = 0. For 
example, Eq. l !118t gives the useful relation 

X m Z" = q mn Z"X m , (m,n) 6 N 2 (122) 

The fundamental relationship between the matrices X and Z is emphasized by the following 
proposition. 

Proposition 6. The unitary matrices X and Z satisfy the q-commutation relation 

[X, Z] q = XL - qZX = (123) 

and the cyclicity relations 

X d = Z d = I d (124) 

In addition, they are connected through 

(F 00 ) + XFoo = Z (125) 
that indicates that X and Z are related by an ordinary DFT transform. 

According to Proposition 6, the matrices X and Z co nstitute a We yl pair (X, Z). Weyl pairs were 
introduced at the beginning of qua ntum mechanics <Weyllll93ll) and used for building opera- 
tor unitary bases jSchwingerl.ll960h . We shall emphasis their interest for quantum information 
and quantum computing in section 4. 

Let x and z be the linear operators associated with X and Z, respectively. They are given by 
x = vqq, z = (v 00 ) f v i =4> xz = v i (126) 



as functions of the operator v ra . Each of the relations involving X and Z can be transcribed in 
terms of x and z. 

The properties of x follow from those of v ra with r = a = 0. The unitary operator x is a shift 
operator when acting on \ j,tn) or \ n) (see {34} and {88)) and a phase operator when acting on 
\ja;00) = |0a;0) (see {47) and {89}). More precisely, we have 

x\j,m) = \j,m@\) x\n) = \nel) (127) 

and 

x|0a;0) = q- K \0a;0) (128) 

The unitary operator z satisfies 

z\j,m) = q hm \j,m) O z|n> = q n \n) (129) 

and 

z|fla;0) = ij _1 |fl«i;0), «i = a 1 (130) 

It thus behaves as a phase operator when acting on \j,m) or \n) and a shift operator when 
acting on \aa;0). 

In view of {128) and {129}, the two cyclic operators x and z (cf . x d = z d = I) are isospectral 
operators. They are connected via a discrete Fourier transform operator (see Eq. <125t ). 
Let us now define the operators 

u ab = x a z b , a,b = 0,l,...,d-l (131) 

The d 2 operators u ab are unitary and satisfy the following trace relation 

tr ((w fl fe) + « fl <(,<) =d& a ,a<h,V ( 132 ) 

where the trace is taken on the d-dimensional space e{d) = e(2j + 1). This trace relation shows 
that the d 2 operators u a f, are pairwise orthogonal operators so that they can serve as a basis 
for developing any operator acting on the Hilbert space e(d). Furthermore, the commutator 
and the anticommutator of u a j, and u^y are given by 

Kb, U aV ] = (q- ba ' - <r flfc ') u a „ b „, a" = a® a', b" = b® b' (133) 

and 

{u ab , u a , b ,} = [q- ha ' + q-" b '^ u a „ b „, a" = a®a', b" = b®b' (134) 

Consequently, [u ab , U^y\ = if and only if ab' Q ba' = and {u ab , u a > b i} = if and only if 
ab' © ba 1 = (l/2)d. Therefore, all anticommutators {u ab ,u a i b i} are different from if d is an 
odd integer. From a group-theoretical point of view, we have the following result. 
Proposition 7. The set {u ab = x"z b : a, b = 0, 1, . . . , d — 1} generates a d 2 -dimensional Lie algebra. 
This algebra can be seen to be the Lie algebra of the general linear group GL(d,C). The subset {u ab : 
a, b = 0, 1, . . . ,d — 1} \ {uqo} thus spans the Lie algebra of the special linear group SL(d,C). 



A second group-theoretical aspect connected with the operators u a ^ concerns a finite group, 
the so-called finit e Heisenberg -Weyl group WH(Z / dZ), known as the Pauli group P d in quan- 
tum information iKiblerl.l2008|) . The set \u a b : a,b = 0,1, d — 1} is not closed under multi- 
plication. However, it is possible to extend the latter set in order to have a group as follows. 
Proposition 8. Let us define the operators w a \, c via 

w abc = ( f u bci a, b,c = 0,1, ... ,d — 1 (135) 

Then, the set {w a b c = q"x b z c : a, b,c = 0,1, ... ,d — 1}, endowed with the multiplication of operators, 
is a group of order d 3 isomorphic with the Heisenberg-Weyl group WH(Z/dZ). This group, also 
referred to as the Pauli group P d , is a nonabelian (for d > 2) nilpotent group with nilpotency class 
equal to 3. It is isomorphic with a finite subgroup of the group U(d) for d even or SU(d) for d odd. 
Proposition 8 easily follows from the composition law 

w abc w a'b'c' = w a"b"c"' a" = a Q a' Q cb' , b" = b © b' , c" = c ® c' (136) 

Note that the group commutator of the two elements w a ^ c and w a iy c i of the group WH(Z/dZ) 
is 

w abc w a'yc\ w abcY ll y w a'yc'Y X = w a"00, a " = bc ' © ch ' (137) 

which can be particularized as 

UabUa'b'iUabrHua'b')- 1 = l""' "" 'l (138) 
in terms of the operators W B f,. 

All this is reminiscent of the group SU(2), the generators of which are the well-known Pauli 
matrices. Therefore, the operators u fl j, shall be referred as generalized Pauli operators and 
their matrices as generalized Pauli matrices. This will be considered further in section 4. 

3.1.5 Link with the cyclic group Q 

There exists an interesting connection between the operator v ra and the cyclic group Q (see 
section 2). The following proposition presents another aspect of this connection. 
Proposition 9. Let R be a generator of Q {e.g., a rotation ofln/d around an arbitrary axis). The 
application 

R" X" : n = 0,l,...,d-l (139) 

defines a d-dimensional matrix representation of C d . This representation is the regular representation 
ofC d . 

Thus, the reduction of the representation {X" : n = 0, 1, . . . ,d — 1} contains once and only 
once each (one-dimensional) irreducible representation 

x («) = (l,<f,...,c i ( i - 1 ) a ) l a = 0,l,...,d-l (140) 

of Q. 



3.1.6 Link with the W M algebra 

Let us define the matrix 



T(„ 1( „ 2 ) = <?i" inj Z">X"S ( Ml ,n 2 )6N 2 (141) 

It is convenient to use the abbreviation 

(n lr n 2 )=n => T (ni#Ill) = T„ (142) 

The matrices T n span an infinite-dimensional Lie algebra. This may be precised as follows. 
Proposition 10. The commutator [T m ,T n ] is given by 

[T m ,T„] = -2/ sin (jm x nj T m+n (143) 

where 

m x n = m-yn.2 — tn2 n l> m + n = (m\ + n.\,mi + n 2 ) (144) 
The matrices T m can be thus formally viewed as the generators of the infinite-dimensional Lie algebra 

Woo. 

The proof of (143) is easily obtained by using (122}. This leads to 

T m T n = q-i mxn T m+n (145) 

which implies (1431. Thus, we get of Lie algebra Woo (or sine algebra) investigated in (Fairlie 
et al., ll990h . 

3.2 Quadratic discrete Fourier transform 
3.2.1 Generalities 

We are now prepared for discussing analogs of the transformations (82) and (84) in the lan- 
guage of classical signal theory. 
Definition 5. Let us consider the transformation 

x = {x m eC:m = Q,l,...,d-\} o y = {y n e C : n = 0,1, . . .,d - 1} (146) 

defined by 

d-l d-l 

yn = E ( F «) m „ *m x m = £ (F ra ) nm y„ (147) 

m=0 n=0 

(F ra )„,„ = * ,^-»W2+(d-i)V4+»[»-(d-l>/q n , m = ,l,...,d-l (148) 
V d 

For a 0, the bijective transformation x <-> y cam be thought of as a quadratic DFT. 
In Eq. (1471 , we choose the matrix F ra as the quadratic Fourier matrix instead of the matrix H ra 
because the particular case r = a = corresponds to the ordinary DFT (see also (Atakishiyev 
et al.,|201fjj)). Note that the matrices F ra and H ra are interrelated via 

(Fra)nm = (H„W n' = d - 1 - n (149) 



Therefore, the lines of F ra in the order 0, 1, . . . , d — 1 coincide with those of H ra in the reverse 
order d - 1, d - 2, . . . , 0. 

The analog of the Parseval-Plancherel theorem for the ordinary DFT can be expressed in the 
following way. 

Theorem 2. The quadratic transformations x y and x' •(-> y' associated with the same matrix F ra , 
with reE and a 6 Z/dZ, satisfy the conservation rule 

i-1 d-l 

Y^Wy'n = E *m x 'm ( 150 ) 

n=0 m=0 

where both sums do not depend on r and a. 
3.2.2 Properties of the quadratic DFT matrix 

In order to get familiar with the quadratic DFT defined by 1 11471 , we now examine some of the 
properties of the quadratic DFT matrix F ra . 

Proposition 11. For d arbitrary, the matrix elements o/F ra satisfies the useful symmetry properties 
( F ra) d -i„ = q (*-V(r+«)n-« e -in(d-l> (F ra ) 0ff , « = 0, 1 d - 1 (151) 

( F ra)„-i a = < ? ( rf - 1 )(''+ fl )/ 2 -«+»« (F ra ) M , n = 1,2, . . .,d - 1, a = 0, 1, . . . ,d - 1 (152) 
zy/tt'c/i can be reduced to the sole symmetry relation 

(Foa)„ei« = '? ( ^ 1)fl/2 ^ + "' , ( F oa), ia , »,« = 0,l,...,d-l (153) 

a;/;en r = 0. 

Proposition 12. For d arbitrary, the matrix F ra is unitary. 

The latter result can be checked from a straightforward calculation. It also follows in a simple 
way from 

(j0L;ra\ip;sb) = (act;r\bp;s) = ((F ra ) f F sb ) a/5 (154) 
It is sufficient to put s = r and b = a in l !154t and to use {48). 

For d arbitrary in addition to be unitary the matrix F ra is such that the modulus of each of its 
matrix elements is equal to 1 / \fd. Thus, F ra can be considered as a generalized Hadamard 
matrix (we adopt here the normalizati on of Hadam ard matrices generally used in quantum 
information and quantum computing jKiblerll2009h ). In the case where d is a prime number, 
we shall prove in section 4 from {154} that the matrix (F ra ) F r b is another Hadamard matrix 
for b a. Similar results hold for the matrix H ra . 
Proposition 13. For d arbitrary, the matrix F ra can be factorized as 

F ra = D ra F, F = F 00 (155) 

where D ra is the d x d diagonal matrix with the matrix elements 

(D ra ) m „ = q Md-m)a/2+(d-lfr/i- m (d-l)r/2 6 ^ n (15fi) 

and F is the well-known ordinary DFT matrix. 

For fixed d, there is one d-multiple infinity of Gaussian matrices D ra (and thus F ra ) distin- 
guished by a 6 Z/dZ and r 6 R. The matrix F was the object of a great number of studies. 



The main properties of the ordinary DFT matrix F are summarized in jAtakishiyev et al],l2010|) . 
Let us simply recall here the fundamental property 

F 4 = I d (157) 

of interest for obtaining the eigenvalues and eigenvectors of F. 
Proposition 14. The determinant o/F ra reads 

detF ra = e ^(rf 2 -i)«/6 detF (158 ) 

where the value o/detF is well-known (Atakishh iev et al\\201(k\Mehta\,\l98~ft) . 
Proposition 15. The trace o/F ra reads 

trF ra = e m ^-^ 2/ ^-^S{u,v,iv) (159) 
v d 

where S(u,v,w) is 

\w\-\ 

S{u,v,w)= £ e m{ukl+vk V w (160) 
Jt=o 

with 



u = 2 — a, v = d(a — r) +r, w = d 



(161) 



(note that v is not necessarily an integer). 

Let us recall that the sum defined by l !160t is a generalized quadratic Gauss sum. It can be 
calculated easily in the situation where u, v and w are integer s such that u and w are mutually 
prime, uw is not zero, and uw + v is even iBerndt et aUll99§) . 

Note that the case a = 2 deserves a special attention. In this case, the quadratic character of 
tr F ra disappears. In addition, if r = we get 



trFn 



Vd 



(162) 



as can be seen from a direct calculation. 

Example 5: In order to illustrate the preceding properties, let us consider the matrix 



?02 



1 



f 1 



\q 5 



1 



<? 
1 
1 



1/ 



(163) 



corresponding to d = 6 (= 



Jn/3 



), r = and a = 2. It is a simple matter of trivial 



calculation to check that the properties given above for F ra are satisfied by the matrix Fq2- 



4. Application to quantum information 



4.1 Computational basis and standard St/ (2) basis 

In quantum information science, we use qubits which are indeed state vectors in the Hilbert 
space C 2 . The more general qubit 

\f 2 ) = c |0) +c x \l), c 6C, cjeC, |c | 2 + | Cl | 2 = l (164) 

is a linear combination of the vectors |0) and |1) which constitute an orthonormal basis 

B 2 = {|0),|1)} (165) 

of C 2 . These two vectors can be considered as the basis vectors for the fundamental irreducible 
representation class of SU(2), in the SU(2) D (1(1) scheme, corresponding to / = 1/2 with 

|0) = |l/2,l/2), |1) = |l/2,-l/2) (166) 

More generally, in d dimensions we use qudits of the form 

d-l d-1 

\ip d )=£ l c„\n), c„eC, n = 0, 1,. . .,d - 1, £|c„| 2 = l (167) 

n=0 n=0 

where the vectors 0), 1 1), . . . , \d — 1) span an orthonormal basis of £r with 

{n\n')=6 lhnl (168) 

By introducing 

j=\{d-\), m = n- X -(A-\), \j,m) = \d - 1 - n) (169) 

(a change of notations equivalent to {68}), the qudits \n) can be viewed as the basis vectors 
\j,m) for the irreducible representation class associated with / of SU(2) in the SU(2) D Lf(l) 
scheme. More precisely, the correspondence between angular momentum states and qudits is 

|0> = |/,;>, |l> = |;,;-l>, \d- 1) = (170) 

where \j,f), \ — 1), ... , \ j, —j) are common eigenvectors of angular momentum operators j 2 
and j z . In other words, the basis (see ll83t), known in quantum information as the compu- 
tational basis, may be identified to the SU(2) D IT(1) standard basis or angular momentum 
basis B 2 j + i (see f3D). We shall see in section 4.2 that such an identification is very useful 
when d is a prime number and does not seem to be very interesting when d is not a prime 
integer. Note that the qudits |0), 1), . . . , \d — 1) are often represented by the column vectors 
(pQ, <p\,..., <pd-i (given by (94)), respectively. 



4.2 Mutually unbiased bases 

The basis B ra given by J701 can serve as another basis for qudits. For arbitrary d, the couple 
(B ra , B d ) exhibits an interesting property. For fixed d, r and a, Eq. dTTb gives 

1 

Vn 6 Z/dZ, Va € Z/dZ : |(n|a*;r)| = -p (171) 

v d 

Equation l|1711 shows that B ra and B ( ; are two unbiased bases. 

Other examples of unbiased bases can be obtained for d = 2 and 3. We easily verify that the 
bases B,o and B rl for d = 2 given by ([58) are unbiased. Similarly, the bases B r Q, B r \ and B r 2 for 
d = 3 given by i61\ are mutually unbiased. Therefore, by combining these particular results 
with the general result implied by \17\\ we end up with three MUBs for d = 2 and four MUBs 
for d = 3, in agreement with N MUB = d + 1 when d is a prime number. The results for d = 2 
and 3 ca n be generalized in the case where d is a prime number. T his leads to the following 
theorem jAlbouy & Kiblerl 120071: lKjbii3. 12008L [Kibler & PlanaAl2006|) . 

Theorem 3. For d = p, with p a prime number, the bases B r Q, B rl , . . ., B ) . p _ 1 , Bp corresponding 
to a fixed value of r form a complete set of p + 1 MUBs. The p 2 vectors \aa;r) or <p(au.;r), with 
a, a = 0,1, ... ,p — 1, of the bases B r Q, B r \, . . . , B rp _i are given by a single formula, namely, Eq. {72j 
or i l93l . The index r makes it possible to distinguish different complete sets ofp + 1 MUBs. 
The proof is as follows. First, according to \17\\ , the computational basis Bp is unbiased with 
any of the p bases B,.q, B r \, ... , B rp _i. Second, we get 

(aa;r\bB;r) = I (172 ) 

or 

(aa;r\bB;r) = - ^ e in{(a-b)B+[p(b-a)+2^-cc)]k}/P (173) 
P Jc=0 



The right-hand side of <173t can be expressed in terms of a generalized quadratic Gauss sum. 
This leads to 

{atx;r\bB;r) = ^S(u,v,w) (174) 



where the Gauss sum S(u, v, w) is given by <160t with the parameters 

u = a — b, v = — (a — b)p — 2(a — B), w = p (175) 

which ensure that uw + v is even. The ge neralized Gauss sum S(w, v, w) in <174t - fl75t can be 
calculated from the methods described in jBerndt et al.l,ll99§) . We thus obtain 

\{aoL;r\bB;r)\ = ^— (176) 
VP 

for all a, b, a, and B in Z/ pZ with b ^ a. This completes the proof. 

Theorem 3 renders feasible to derive in one step the (p + l)p qupits (i.e., qudits with d = p a 
prime integer) of a complete set of p + 1 MUBs in C p . The single formula $72\ or (93) , giving 
the p 2 vectors \a&; r) or tp (ax; r), with a, a = 0,1, ... ,p — 1, of the bases B,-o, B rl , . . . , B rp _ 1 , is 
easily codable on a classical computer. 



Example 6: The p = 2 case. For r = 0, the p + 1 = 3 MUBs are 



d . |0) + |1> |0) 
£>oo ■ 



V2 ' y/l 

. |0> -i|l> . |o)+i|i) n _ 

B 2 : |0), |1) 

cf. (59), (60), jlOOt and jlOlt . The global factors —1 in Bqo and ±z in B i arise from the general 
formula {72); they are irrelevant for quantum information and can be omitted. 
Example 7: The p = 3 case. For r = 0, the p + 1 = 4 MUBs are 



300 



=01 



302 



|0) + |1) + |2) 


< ? 2 |0)+< ? |1) + |2) 


q\0)+q 2 \l} + \2) 




V3 


V3 


q\0) +q\\) + \2) 


|0)+q 2 |l) + |2) 


q 2 \0) + \l) + \2) 


V3 


V3 


V3 


q 2 \0)+q 2 \l) + 


2) q\0) + |1) + |2> 


|0>+<7|1> + |2> 


V3 


V3 


V3 



(178) 



B 3 : |0), |1), \2) 
with q = e i27z/3 , cf. (62), (63), (M), (106) . (108) and (TlOt . 

As a simple consequence of Theorem 3, we get the following corollary which can be derived 
by combining Theorem 3 with Eq. (154). 

Corollary 2. For d = p, with p a prime number, the p x p matrix (F ra ) + F r i, with b ^= a (a,b = 
0, 1, . . . , p — 1) is a generalized Hadamard matrix. 

Going back to arbitrary d, it is to be noted that for a fixed value of r, the A + 1 bases B r Q, B rl , . . ., 
Brd-1/ Srf do not provide in general a complete set of d + 1 MUBs even in t he case wher e d is 
a power p e with e > 2 of a prime integer p. However, it is possible to show (Kiblerl.l2009|) that 
the bases Bo„, B 0a(S i and B^ are three MUBs in C d , in agreement with N MUB > 3. Therefore for 
d arbitrary, given two Hadamard matrices F ra and F s j,, the product F ra + F s i, is not in general a 
Hadamard matrix. 

In the case where d is a power p e with e > 2 of a prime integer p, tensor products of the 
unbiased bases B,.q, B r \, . . ., B rp _i can be used for generating p e + 1 MUBs in dimension 
d = p c '. This can be illustrated with the following example. 

Example 8: The d = 2 2 case. This case corresponds to a spin j = 3/2. The application of 
(45) or (72) yields four bases Bq q (a = 0, 1, 2, 3). As a point of fact, the five bases Bqo, Bqi, Bq2, 
B03 and B 4 do not form a complete set of d + 1 = 5 MUBs (d = 4 is not a prime number). 
Nevertheless, it is possible to find five MUBs because d = 2 2 is the power of a prime number. 
This can be achieved by replacing the space e(4) spanned by 

B 4 = {|3/2,m) : m = 3/2,1/2,-1/2,-3/2} or {\n) : n = 0,1,2,3} (179) 

by the tensor product space e(2) ® e(2) spanned by the canonical or computational basis 

B 2 ®B 2 = {|0) ® |0), |0) ® |1), |1) ® |0), |1) ® |1)} (180) 



The space e(2) ® e(2) is associated with the coupling of two spin angular momenta j\ = 1/2 
and ]2 = 1/2 or two qubits (in the vector u®v, u and v correspond to j\ and ]i, respectively). 
Four of the five MUBs for d = 4 can be constructed from the direct products 

\ab : /up) = \aa;Q) ® |fo/3;0) (181) 

which are eigenvectors of the operators 

Wab = »0fl ® »0f> (182) 

(the operators vq„ and i; ij refer to the two spaces e(2), the vectors of type |aa;0) and &B;0) 
are given by the master formula J72b for d = 2). Obviously, the set 

Bo«ob = {l«&:<*j6) :«,/3 = 0,l} (183) 

is an orthonormal basis in C 4 . It is evident that Bqooo an d Bqioi are two unbiased bases since 
the modulus of the inner product of 1 00 : a/3) by 1 11 : ct'j6') is 

(00 : ajS|ll : a'j6')| = |(0a;0|ln;';0)(0/3;0|l^;0)| = ^=^1= = 4= (184) 

v2 v2 v4 

A similar result holds for the two bases Bqooi and Bqioq. However, the four bases Boooo B0101, 
Bqooi and Bqioo are not mutually unbiased. A possible way to overcome this no-go result is to 
keep the bases Bqooo an d Bqioi intact and to re-organize the vectors inside the bases Bqooi and 
Bqioo m order to obtain four MUBs. We are thus left with four bases 

W 00 = Boooo, Vf u = Bqioi, W i, W w (185) 
which together with the computational basis B4 give five MUBs. Specifically, we have 

W 00 = { 00 : aj8) : a,6 = 0,1} (186) 

Wn = {|11 : ctp) : a,fi = 0,1} (187) 

W 01 = {A|01 : tip) +}i\0\ : a©16©l) : a,6 = 0,1} (188) 

W 10 = {A|10:a6)+/(|10:a©16©l) :a,6 = 0,l} (189) 



where 



A=l^, = (190) 



As a resume, only two formulas are necessary for obtaining the d 2 = 16 vectors for the bases 
W ab , namely 

W 00 ,Wn : \aa:ap) (191) 
W i,W 10 : A|flflffil : aB) +^|flflffil : a© 16© 1) (192) 

for all a, a and B in Z / 2Z. The five MUBs are listed below as state vectors and column vectors 
with 

|0) = |; = 1/2, m = 1/2) or (l j , |1) = |j = 1/2, m = -1/2) or (J J (193) 



The canonical basis: 



or in column vectors 



10)810), |0)®|1), |1)®|0), |1)®|1) 






















1 






















1 







W 




W 











The Wqo basis: 

1 00 : 00) = + 

1 00 : 10) = - 
or in developed form 



0) + |1) |0) + |1) 



o> - |i> „ |Q) + |i) 



1 00 : 01) = 



1 00 : 11) = + 



10) + 11) 10) -11) 

|0)-|1) n |0)-|1) 
V2 



V2 



1 00 : 00) = +-(|0) ® |0) + |0) ® |1) + |1) ® |0) + |1) ® |1)) 
1 00 : 01) = --(|0) <g> |0) - |0) ® |1) + |1) ® |0) - |1) ® |1)) 
1 00 : 10) 



(|0) ® |0) + |0) ® |1) - |1) ® |0) - |1) ® |1)) 



1 00 : 11) = +-(|0) ® |0) - |0) ® |1) - |1) ® |0) + |1) ® |1)) 



or in column vectors 

1 
2 

The Wii basis: 



/A 
1 
1 

w 



1 

v-v 



-1 

v-v 



-1 
V i / 



|11 : 00) = - 



|0)-f|l) ^ |0> — z|X> 



|ll:10)=+M+M 
n/2 



|Q)-'|i) 



|11 : 01) = + 



|ll:ll ) = -'°> + ;' 1 > 



\0)-i\l) CT |0) + »|1) 
|0)+i|l) 



or in developed form 

|11 : 00) = -i(|0)®|0) - i|0)g> |l)-i|l)®|0) - |1) ® |1)) 
11 : 01) = +i(|0)®|0) +i|0)® |l)-z|l)®|0) + |1) ® |1)) 
|11 : 10) = +i(|0)®|0) -i|0)® |l)+z|l)®|0) + |1) ® |1)) 
|11 : 11) = -i(|0)® |0) +i|0)® |l)+i|l)® |0) - |1) ® |1)) 



(194) 



(195) 



or in column vectors 





f 1 ^ 




(1\ 
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1 


—i 


1 


i 


1 


—i 


1 


2 


—i 


' 2 


—i 


' 2 


i 


2 


/ 



















(196) 



The Wni basis: 



A|01 : 00) + p|01 : 11) = +f( 



p|01 : 00) +A|01 : 11) 



A|01 : 01) +}i\01 : 10) 



fi\0l : 01) +A|01 : 10) = +A 



|0) + |1) |0)-/|1) . |0> — |X> |0)+z|l) 
A- 



V2 



V2 



M+ji> 8 M-i!i> 

y/2 V2 



V2 

Q)-|i) 

V2 



|0) + |1) ^ |0)+i|l) , J0>-|1> 



V2 



V2 



|0) + |1) ^ |0)+i|l) 



V2 



V2 

Q)-|i) 

V2 



|0) + '|1) 

|Q)-'|i) 

|0)-'|1) 
V2 



or in developed form 



A|01 : 00) + p|01 : 11) = +^(10) ® |0) - |0) ® |1) - ® |0) - i 



1>®|1» 
1>®|1» 

(|0) ® |0) + |0) ® |1) - ® |0) + t|l) ® |1)) 

1)0 |1» 



;/ 101 : 00) +A|01 : 11) = +-(|0) ® |0) + |0) ® ® |0) -t 

A|01 : 01) +/(|01 : 10) = 



ii|01 : 01) +A|01 : 10) = -^(10) ® |0) - |0) ® |1) + ® |0) + i 



or in column vectors 



The W w basis: 
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-1 


— i 


' 2 
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' ~2 




' ~2 
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\ * J 
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A|10 


00) + 1(|10 


11) 


- i 
V2 




|0) + |1) 


A |0)+i 


1) 


jo) -|i) 


n\w 


00) +A|10 


11) 


JO)-i 




|o) + |i) 


|0)+i 


— e 


, |0)-|1> 
V2 


A|10 


01) +ji\W 


10) 


10) -i 




|0> — |X> 


+a'°> + ; 




, |0> + |1> 


F|!0 


01) +A|10 


10) 


= + a'°>-J 


*>« 


, |0> — |X> 


\Q) + i 
~ }l V2 


1) 


a |o> + |i> 
V2 



or in developed form 



A|10 


00) +}i 


10 


11) 


= +2(|0>6 


3 10) - i 


0)g 


Hi) 


-|l)c 


5 |0) - 






jt|10 


00) + A 


10 


11) 


= +j(|0>6 


3 10) +i 


0)c* 


)|i> 


+ |i> c 


3 |0) - 


|1> eg 




A|10 


01)+/' 


10 


10) 


= ~(|o>e 


3 10) +i 


0)c* 


5|1> 


-|1>C 


3 |0) + 


|i)e 


9|1» 


p\W 


01) + A 


10 


10) 


= ~(|o>e 


3 \0)-i 


0)c* 


)|i> 


+ |i> c 


3 |0) + 


|1> eg 





or in column vectors 



/ 1 \ 

-i 

-1 



/1\ 

i 

1 



/ 1 ^ 

i 


i 


—i 


-1 


' ~2 


1 


V i J 







(198) 



The five preceding bases are of central importance in quantum information for expressing any 
qu quart or quartic (corresponding to d = 4) in terms of qudits (corresponding to d = 2). It 
is to be noted that the vectors of the Wqo and W-yi bases are not intricated (i.e., each vector is 
the direct product of two vectors) while the vectors of the W 01 and W 10 bases are intricated 
(i.e., each vector is not the direct product of two vectors). To be more precise, the degree of 
intrication of the state vectors for the bases Wqo, W 11; W 01 and W w can be determined in the 
following way. In arbitrary dimension d, let 



If) = EE%W 

/c=0 /=0 



\l) 



(199) 



be a double qudit state vector. Then, it can be shown that the determinant of the d x d matrix 
A = (a k i) satisfies 



< I detAI < 



d d 



(200) 



as discussed in jAlbouyl.l2009T) . The case det A = corresponds to the absence of global intri- 
cation while the case 



det A 



(201) 



corresponds to a maximal intrication. As an illustration, we obtain that all the state vectors for 
Wqo and Wi\ are not intricated and that all the state vectors for Wpj and Wig are maximally 
intricated. 

Generalization of dl9H and 1 11921 can be obtained in more complicated situations (two qupits, 
three qubits, . . . ). The generalization of 1 11911 is immediate. The generalization of J1921 can be 
achieved by taking linear combinations of vectors such that each linear combination is made 
of vectors corresponding to the same eigenvalue of the relevant tensor product of operators 
of type v 0a . 



4.3 Mutually unbiased bases and Lie agebras 
4.3.1 Generalized Pauli matrices 

We now examine the interest for quantum information of the Weyl pair (X, Z) introduced in 
section 3.1.4. The linear operators corresponding to the matrices X and Z are known in quan- 
tum information and quantum computing as shift and clock operators, respectively. (Note 
however that for each of the operators x and z, the shift or clock character depends on which 
state the operator acts. The qualification adopted in quantum information and quantum com- 
puting corresponds to the action of x and z on the computational basis Bj.) For d arbitrary, 
they are at the root of the Pauli group P^, a finite subgroup of U(d) (see section 3.1.4). The 
normaliser of in U(d) is a Clifford-type group in d dimensions noted CZ d . More precisely, 
CZd is the set {U 6 !I(d)|UP ( jU + = Pd} endowed with matrix multiplication (the elements of 
Ptf being expressed in terms of the matrices X and Z). The Pauli group P$, as well as any other 
invariant subgroup of C/^, is of considerable importa nce for describing quantum errors and 
quantum fault toleran ce in quantum computing (see jHavhcek & SanigaTl2008l;|Planaj,l201Ct 
IPlanat & Kiblerl l2010|) and references therein for recent geometrical approaches to the Pauli 
group). These concepts are very important in the case of n-qubit systems (corresponding to 
d = 2"). 

The Weyl pair (X, Z) turns out to be an integrity basis for generating the set {X"Z b : a,b 6 
Z/dZi of d 2 generalized Pauli matrices in d dimensions (see for instance (Bandyopadhyay et 
al., l2002tlGottesman et al.l,l200ll;lKiblerll2008tlLawrence et aUl2002t|Pittenger & Rubinl,l2004 
in th e context of MUBs and dBalian & Itzvksonl,ll986HPatera & Zassenhauslll988l: St'ovlcek & 



Tolar,ll98j in group-theoretical contexts). As seen in section 3.1.4, the latter set constitutes 



a basis for the Lie algebra of the linear group GL(d,C) (or its unitary restriction U(d)) with 
respect to the commutator law. Let us give two examples of these important generalized Pauli 
matrices. 

Example 9: The d = 2 case. For d = 2 / = 1/2 (=> q = —1), the matrices of the four 
operators u a y with a, b = 0, 1 are 

I 2 = X°Z°, X = X 1 Z°, Y = X l Z l , Z = X°Z 1 (202) 

or in terms of the matrices Vg a 

I 2 = (V 00 ) 2 , X = V 00 , Y = Voi, Z = (V 00 ) + V 01 (203) 

In detail, we get 

*-G 9- x =(? ;)• Y =(? "o 1 )' z =G -0 

Alternatively, we have 

I 2 = a , X = a x , Y = -icy, Z = a z (205) 
in terms of the usual (Hermitian and unitary) Pauli matrices c , a~ x , a,j and a z 

—(J 0- »-(! i)- ■*-(! o'). -ft -0 <-» 



The approach developed here leads to generalized Pauli matrices in dimension 2 that differ 
from the usual Pauli matrices. This is the price one has to pay in order to get a systematic 



generalization of Pauli matrices in arbitrary dimension. It should be observed that the com- 
mutation and anti-commutation relations given by <133t and dl34j with d = 2 correspond 
to the well-known commutation and anti-commutation relations for the usual Pauli matrices 
transcribed in the normalization X 1 Z° = a x , X 1 Z 1 = —icy, X°Z 1 = cr z . 

From a group-theoretical point of view, the matrices I2, X, Y and Z can be considered as gener- 
ators of the group IT (2) . On the other hand, the Pauli group P2 contains eight elements; due to 
the factor — i in Y = — ia y , the group P2 is isomorphic to the group of hyperbolic quaternions 
rather than to the group of ordinary quaternions. 

In terms of column vectors, the vectors of the bases Bqo, Box an d B2 (see i!77\ ) are eigenvevtors 
of a x , ay and a z , respectively (for each matrix the eigenvalues are 1 and —1). 
Example 10: The d = 3 case. For d = 3 -O- j = 1 (=> q = e l2n ^), the matrices of the nine 
operators u a f, with a, b = 0, 1, 2, viz., 

X°Z° = I 3 , X^ = X, X 2 Z° = X 2 

X°Z 1 = Z, X°Z 2 = Z 2 , X 1 Z 1 = XZ (207) 
X 2 Z 2 , X 2 Z : = X 2 Z, X 1 Z 2 = XZ 2 

are 

/1 o\ /o 1 0\ 

x 2 = 






/l 0' 
Z = q 
\0 q 2 , 



a o^ 
q 2 

^0 qj 



XZ 



'0 q 0^ 
q 2 
a 0> 



(208) 



1 , XZ 2 = [ 





2 



0; 



The generalized Pauli matrices J208t differ from the Gell-Mann matrices used in elementary 
particle physics. They constitute another extension of the Pauli matrices in dimension d = 3 
of interest for the Lie group !i(3) and the Pauli group P3. 

In terms of column vectors, the vectors of the bases Bqo, Bqi, B02 and B3 (see J178t ) are eigen- 
vectors of X, XZ, XZ 2 and Z, respectively (for each matrix the eigenvalues are 1, q and q 2 ). 



4.3.2 MUBs and the special linear group 

In the case where d is a prime integer or a power of a prime integer, it is known that the set 
{X"Z b : a, b = 0, 1, . . . ,d — 1} of cardi nality d 2 can be partitioned in to d + 1 subsets containing 
each d — 1 commuting matrices (cf . jBandyopadhyay "et mnooi). Let us give an example 
before going to the case where d is an arbitrary prime number. 

Example 11: The d = 5 case. For d = 5, we have the six following sets of four commuting 
matrices 



V = {01,02,03,04}, Vi = {10,20,30,40} 
V 2 = {11,22,33,44}, V 3 = {12,24,31,43} 
V 4 = {13,21,34,42}, V 5 = {14,23,32,41} 



(209) 



where ab is used as an abbreviation of X a Z b . 

Proposition 16. For d = p with p a prime integer, the p + 1 sets of p — 1 commuting matrices are 
easily seen to be 



V 


= {X°Z" : 


a = 1, 2, . . . , p - 


1} 






Vi 


= {X«Z° : 


a = 1, 2, . . . , p - 


1} 






v 2 


= {X fl Z fl : 


a = l,2,...,p- 


1} 






v 3 


= {X fl Z 2fl 


: a = 1, 2, . . . , p - 


-1} 
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p-1 


= {x fl z(p 


- 2 >" :a = l,2,... 


•v- 


1} 




V,, 


= {x a z(p- 


: a = 1,2,... 


>v- 


1} 





Each of the p + 1 sets Vq, Vi, . . . , V» caw be put fw a one-to-one correspondance with one basis of the 
complete set ofp + 1 MUBs. In fact, Vq is associated with the computational basis B p ; furthermore, in 
view of 

Vo a eV„ +1 , 8 = 0,1,. ..,p-l (211) 

it follows that V\, Vz, ■ ■ ., V p are associated with the p remaining MUBs Bqq, Bq\, . . ., Bo p _i, respec- 
tively. 

Keeping into account the fact that the set {X a Z b : a, b = 0, 1, . . . ,p - 1} \ {X°Z } spans the Lie 
algebra of the special linear group SL(p,C), we have the next theorem. 

Theorem 4. For d = p with p a prime integer, the Lie algebra sl(p,C) of the group SL(p,C) can be 
decomposed into a sum (vector space sum indicated by l±l) of p + 1 abelian subalgebras each of dimension 
p — 1, i.e., 

s/(p,C) ~ » ttl z»x tbl . . . ttJ Op (212) 

where the p + 1 subalgebras vq, V\, . . .,v p are Cartan subalgebras generated respectively by the sets 
Vq, Vi, . . . , Vp containing each p — 1 commuting matrices. 

The latter result can be extended when d = p e with p a prime integer and e an integer (e > 2): 
ther e exists a decomposition of sl(p c ,C) into p e + 1 abelian sub algebras of dimension p e — 1 
(cf. jBoykin et al.l.l2007l:lkiblerl.l2009l;rPatera & Zassenhaud.ll98£k 



5. Conclusion 

The quadratic discrete Fourier transform studied in this chapter can be considered as a two- 
parameter extension, with a quadratic term, of the usual discrete Fourier transform. In the 
case where the two parameters are taken to be equal to zero, the quadratic discrete Fourier 
transform is nothing but the usual discrete Fourier transform. The quantum quadratic discrete 
Fourier transform plays an important role in the field of quantum information. In particular, 
such a transformation in prime dimension can be used for obtaining a complete set of mutu- 
ally unbiased bases. It is to be mentioned that the quantum quadratic discrete Fourier trans- 
form also arises in the determination of phase operators for the groups S U(2) and SU(1, 1) 
in connection with the repre sentations of a generalized oscillator algebra jAtakishiyev et all 
l2010l: iDaoud & Kiblerl.l2O10|) . As an open question, it should be worth investigating the re- 
lation between the quadratic discrete Fourier transform and the Fourier transform on a finite 
ring or a finite field. 
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Appendix: Wigner-Racah algebra of SU(2) in the {j 2 ,x} scheme 

In this self-contained Appendix, the bar does not indicate complex conjugation. Here, com- 
plex conjugation is denoted with a star. 

The Wigner-Racah algebra of the group SU(2) in the SU(2) D or {j 2 ,jz} scheme is 

well known. It corresponds to the use of bases of type B 2 y + i resulting from the simultaneous 
diagonalization of the Casimir operator j 2 and of the Cartan generator j z of SU(2). Any change 
of basis of type 



\hv) = E \i> m ){iMi>F) 



(213) 



(where for fixed ; the elements (j,m\j,fi) define a (2; + 1) x (2; + 1) unitary matrix) leads 
to another acceptable scheme for the Wigner-Racah algebra of SU(2). In this scheme, the 
matrices of the irreducible representation classes of SU(2) take a new form as well as the 
coupling coefficients (and the associated 3 — jm symbols). For instance, the Clebsch-Gordan 
or coupling coefficients (/i/2 ra i' M 2 \jm) ar e simply replaced by 



h h i 

{hhnMil 1 ) = E E E (hh m i m 2\jm) 

mi=-;'i m z =-j z m=-j 

(h> m i\h>Hi)* {k> m 2\)2>V2)* {iMi>v) 



(214) 



when passing from the {jm} quantization to the {;/(} quantization while the recoupling co- 
efficients, and the corresponding 3(n — 1) — ; symbols, for the coupling of n (n > 3) angu- 
lar momenta remain invariant. The adaptation to the {jji} quantization scheme afforded by 
Eq. {213} is transferable to SU(2) irreducible tensor operators. This yields the Wigner-Eckart 
theorem in the {j}i} scheme. 

We give here the basic ingredients for developing the Wigner-Racah algebra of SU (2) in the 
{j Z > v 0o} or {j 2 r x } scheme. For such a scheme, the vector is of the form |/a;00) so that 
the label /( can be identified with a. Thus, the inter-basis expansion coefficients (]',m\j, ft) are 



(i,m\ioi; 00) 



Jj+m)a 



exp 



2 m 



2/ + 



-(j + m)a 



(215) 



with m = ),) 
formation d45 



- 1,...,- 
with r 

3 oo- 



-j and a = 0, 1, ... , 2j. Equation J215t corresponds to the unitary trans- 
= a = 0, that allows to pass from the standard basis B 2 y + j to the non- 
standard basis Bnn- Then, the Clebsch-Gordan coefficients in the {j 2 , i>oo} scheme are 



I h n ft 

(w2«i« 2 |/3*3) = V(2/l+1)(2/2+1)(2/3+1) m E_. im E_. 2m E_. 3 

(9l )-(A+mi)«i {q2 )-(i2+»^ 2 ((73) (;3+m 3 )« 3 ( hhmi m 2 \j 3 m 3 ) 



(216) 



where the various are given in terms of fa by 



qk = exp (wri) ' k = 1 ' 2 ' 3 (217) 

The symmetry properties of the coupling coefficients (71/2*1*2 73 a 3) cannot be expressed in a 
simple way (except the symmetry under the interchange j\X\ <-> 72*2)- Therefore, it is inter- 
esting to introduce the following / symbol through 

y(h h h\ = I f f f 

V a l "2 0C 3 J 7(2 ;i + 1) (2; 2 + 1) (2/ 3 + 1) m tL h mi =-p_ m3 = -, 3 

{qi) -(n+»>^ (?2 )-(ft+*2)«2 (?3) -(A+».)«a (fl & (21g) 



where the 3 — /m symbol on the right-hand side of {218} is an ordinary Wigner symbol for the 
SU(2) group in the {j 2 ,jz} scheme. (The / symbol is to the {j 2 , x} scheme what the V symbol 
of Racahis to the {j 2 ,jz} scheme.) The/symbol exhibits the same symmetry properties under 
permutations of its columns as the 3 — jm Wigner symbol (identical to the V Racah symbol): 
Its value is multiplied by (— \}h+h+fr under an odd permutation and does not change under 
an even permutation. In contrast to the 3 — jm symbol, not all the values of the / symbol are 
real. In this respect, the / symbol behaves under complex conjugation as 

Other properties (e.g., orthogonality properties, connection with the Clebsch-Gordan coeffi- 
cients and the Herring-Wigner tensor, etc.) of the / symbol an d its relations with 3(n — 1) — j 
symbols for n > 3 can be derived along the lines developed in |Kiblerl.ll968|) . 
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